Three-spin interactions in optical lattices and criticality in cluster Hamiltonians 
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We demonstrate that in a triangular configuration of an optical lattice of two atomic species a variety of 
novel spin- 1/2 Hamiltonians can be generated. They include effective three-spin interactions resulting from the 
possibility of atoms tunneling along two different paths. This motivates the study of ground state properties of 
various three-spin Hamiltonians in terms of their two-point and n-point correlations as well as the localizable 
entanglement. We present a Hamiltonian with a finite energy gap above its unique ground state for which the 
localizable entanglement length diverges for a wide interval of applied external fields, while at the same time 
the classical correlation length remains finite. 
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The combination of cold atom technology with optical lat- 
tices [1,2] gives rise to a variety of possibilities for construct- 
ing spin Hamiltonians [3,4]. This is particularly appealing 
as the high degree of isolation from the environment that can 
be achieved in these systems allows for the study of these 
Hamiltonians under idealised laboratory conditions. In paral- 
lel, techniques have been developed for minimising imperfec- 
tions and impurities [5,6] in the implementation of the desired 
structures and for their subsequent probing and measurement 
[7]. These achievements permit the experimental investiga- 
tion of Hamiltonians that are of interest in areas such as quan- 
tum information or condensed matter physics with the added 
advantage of a remarkable freedom in the choice of exter- 
nal parameters. Presently, attention both in condensed matter 
physics and in cold atom research is focusing on two-spin in- 
teractions as these are most readily accessible experimentally. 
However, the unique experimental capability provided by cold 
atom technology allows us to relax this restriction. Here we 
demonstrate that cold atom technology provides a laboratory 
to generate and study higher order effects such as three-spin 
interactions that give rise to unique entanglement properties. 

The present work serves two purposes. Firstly, it demon- 
strates that in a two species Bose-Hubbard model in a triangu- 
lar configuration a wide range of Hamilton operators can be 
generated that include effective three-spin interactions. They 
result from the possibility of atomic tunneling through dif- 
ferent paths from one vertex to the other. This can be ex- 
tended to a one dimensional spin chain with three-spin inter- 
actions. Secondly, we take this novel experimental capabil- 
ity as a motivation to study unique ground state properties 
of Hamiltonians that include three-spin interactions. In this 
context one can study possible quantum phase transitions by 
considering both the classical correlation properties as well 
as the entanglement properties of these systems. Specifically 
we consider the so-called cluster Hamiltonian and its ground 
state, the cluster state which has previously been shown to 
play an important role as a resource in the context of quantum 
computation [20]. Subject to an additional Zeeman term the 



combined Hamiltonian possesses a finite energy gap above its 
unique ground state in a finite parameter range, hence exhibit- 
ing no critical behaviour in the classical correlations in that 
regime. We shall show that at the same time it exhibits a crit- 
ical behaviour in its entanglement properties due to its three 
spin- 1/2 interaction term. This is manifested by a diverging 
entanglement length of the localizable entanglement [8]. Our 
example demonstrates that divergence in entanglement prop- 
erties are not necessarily related to the existence of classical 
critical points, the latter giving a rather incomplete description 
of the long-range quantum correlations against popular belief 
[9]. A related example was arrived at independently in [10]. 

Consider an ensemble of ultracold bosonic atoms confined 
in an optical lattice formed by several standing wave laser 
beams [3,4,11]. Each atom is assumed to have two relevant 
internal states, denoted with the index a — a,b, which are 
trapped by independent standing wave laser beams differing 
in polarisation. We are interested in the regime where the 
atoms are sufficiently cooled and the periodic potential is high 
enough so that the atoms will be confined to the lowest Bloch 
band and the low energy evolution can be described by the two 
species Bose-Hubbard Hamiltonian [11]. The tunneling cou- 
plings and the coUisional couplings Uaa' can be widely 
varied by adjusting the amplitude of the lattice laser fields. 
For the generation of the multi-particle interactions discussed 
here we require large collisional couplings in order to have a 
significant effect within the decoherence time of the system. 
This can be achieved experimentally by Feshbach resonances 
[12-14], for which first theoretical [15] and experimental [16] 
advances are already promising. 

Let us begin by considering the case of only three sites in 
a triangular configuration (see Figure 1) with tunneling cou- 
pling activated between all three of them. We are interested 
in the regime where the tunneling couplings are much smaller 
than the collisional ones, ^ Uaa' which corresponds to 
the Mott insulating phase and we demand that we have on 
average one atom per lattice site. Hence, the basis of states 



of site i can be defined by 1 1) 
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1 1) = |n° = 0,71^ = 1), where n° and are the num- 
ber of atoms in state a or 6 respectively. It is possible to ex- 
pand the resulting evolution generated by the Bose-Hubbard 
Hamiltonian in terms of the small parameters J"' /Uaa'- In 
an interaction picture with respect to the collisional Hamil- 
tonian, i?^") = i X^mo-' ^<^c^'"lfT«lfT'°^»o-'««c^' obtains 
the effective evolution from the perturbation expansion up 
to the third order with respect to the tunneling interaction, 

V = - E^aiJ>l<^^+l'r + H-C), givCn by 
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The indices a, (3 refer to states with one atom per site while 7, 
5 refer to states with two or more atomic populations per site, 
are the eigenvalues of the collisional part, H^^\ while we 
neglected fast rotating terms effective for long time intervals 
[17]. Written explicitly in terms of spin operators we obtain 
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+A(3)afa^^a^, + AW(afa^ia?V2 + ^t<Tt+i^^S\ ■ (2) 
The couplings A''-' are given as an expansions in J'^ /Uaa' by 
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The local field B can be arbitrarily tuned by applying appro- 
priately detuned laser fields while we need to compensate for 
single particle phase rotations of the form '^i with 

J*^^ , 9 J*^ J'' , J^^ , 9 J'' J'' , 
B^ = (2 H \ ) H (2 H h ). 

Uaa "^Uaa Uab Ubb 2Ubb Uab 

One can isolate different parts from Eq. (2), each one in- 
cluding a three-spin interaction term, by varying the tunneling 
and/or the collisional couplings appropriately so that particu- 
lar A*-*' terms such as the two spin interactions vanish, while 
others can be varied freely. 

By employing additional Raman transitions in such a way 
as to couple the states a and b during tunneling it is possi- 
ble to obtain variations of the above Hamiltonian [4]. In- 
deed, Raman transitions can activate tunneling of the states 



1+) = ^(|T) + U)). while the tunnehng of the states 
|— ) = -^(11) — 14)) is obstructed. Hence, it is possible to gen- 
erate different coefficients in front of the two spin interaction 
terms crfaf^i or erf ct^'^j^ as they are the diagonal matrices in 
their corresponding basis |+) or |— ). Considering its effect on 
the three-spin interaction it is possible to generate additional 
terms of the form crf^^crf^j °r '^i with coupHngs 

similar to A'-'^-'. Note that the effective spin interactions pro- 
duced by Raman transitions do not preserve the number of 
particles in each of the species. 

In particular, we are interested in obtaining a whole 
chain of triangles in a zig-zag one dimensional pattern as 
in Fig. 1. Indeed, with this configuration we can extend 
from a single triangle to a whole triangular ladder. Nev- 
ertheless, a careful consideration of the two spin interac- 
tions shows that terms of the form afafj^2 ^Iso appear, due 
to the triangular configuration (see Fig. 1). Hamiltoni- 
ans involving nearest and next-to-nearest neighbour interac- 
tions are of interest in their own right (see e.g. Chapter 
14 of [9] and [18]), but we will not address these systems 
here. It is possible to introduce a longitudinal optical lat- 
tice with half of the initial wave length, and an appropri- 
ate amplitude such that it cancels exactly those interactions 
generating finally chains with only neighbouring couplings. 



iH-1 




FIG. 1 . The one dimensional chain constructed out of equilateral 
triangles. Three-spin interaction terms appear, e.g. between sites i, 
i + 1 and i + 2 as, for example, tunneling between i and i + 2 can 
happen through two different paths, directly and through site i + 1. 
The latter resulting into an interaction between i and i + 2 that is 
controlled by the state of site i + 1. 



In a similar fashion it is possible to avoid generation of the 
term erf (Tf_^2 + cr|'cr|'^2 by deactivating the longitudinal tun- 
nelling coupling in one of the modes, e.g. the a mode which 
deactivates the corresponding exchange interaction. We are 
particularly interested in three-spin interactions and would 
Uke to isolate the chain term Yl,Mi'^t+i'^i+2 + 
from the A*^"'^ term (see Hamiltonian (2)) that includes in ad- 
dition all the possible triangular permutations. To achieve that 
we could deactivate now the non-longitudinal tunnelling for 
one of the two modes, e.g. the a-mode. With the above pro- 
cedures we can finally obtain a chain Hamiltonian as in (2) 
where the summation runs up to the total number N of the 
sites. A variety of different Hamiltonians could be generated 
by different combinations of the above techniques. 

In the past Hamiltonians describing three-spin interactions 
have been of limited interest [19] as they were difficult to 
implement and control experimentally. The above results 
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demonstrate that Hamiltonians with three-spin interactions 
can be implemented and controlled across a wide parameter 
range. One may suspect that ground states of three-spin inter- 
action Hamiltonians exhibit unique properties as compared to 
ground states generated merely by two-spin interaction. This 
motivates the study of the properties of the ground state of 
a particular three-spin Hamiltonian for different parametric 
regimes. Possible phase transitions induced by varying these 
parameters are explored employing two possible signatures of 
critical behaviour that are quite different in nature. In partic- 
ular, new critical phenomena in three-spin Hamiltonians that 
cannot be detected on the level of classical correlations will 
be demonstrated. 

(i) A traditional approach to criticality of the ground state 
studies two-point correlation functions between spins 1 and L, 
given by C"^ = (cr^cr^) — (cr")((T^), for varying L, where 
a,(3 = X, y, z. These two-point correlations may exhibit two 
types of generic behaviours, namely (a) exponential decay in 
L, i.e. the correlation length ^, defined as 

lim y log Cit, (3) 

is finite or, (b), power-law decay in L, i.e. C"^ ^ L~'^ for 
some q, which implies an infinite correlation length ^ indicat- 
ing a critical point in the system [9]. 

(ii) While the two-point correlation functions C"^ are a 
possible indicator for critical behaviour, they provide an in- 
complete view of the quantum correlations between spins 1 
and L. Indeed they ignore correlations through all the other 
spins by tracing them out. Already the GHZ state \GHZ) = 
(1 000) + I 111)) /a/2 shows that this looses important infor- 
mation. Tracing out particle 2 leaves particles 1 and 3 in an 
unentangled state. However, measuring the second particle in 
the cra;-eigenbasis leaves particles 1 and 3 in a maximally en- 
tangled state. Therefore one may define the localizable entan- 
glement i^l'^'^-' between spins 1 and L as the largest average 
entanglement that can be obtained by performing optimised 
local measurements on all the other spins [8]. In analogy to 
Eq. (3) one can define the entanglement length 

C^'^ lim ylogi;fr'. (4) 

It is an interesting question whether criticality according 
to one of these indicators impUes criticality according to the 
other. The localizable entanglement length is always larger 
than or equal to the two-point correlation length and indeed, 
it has been shown that there are cases where criticaUty be- 
haviour can be revealed only by the diverging localizable en- 
tanglement length while the classical correlation length re- 
mains finite [10]. Such behaviour is also expected to appear 
when we consider particular three-spin interaction Hamiltoni- 
ans. To see this consider the Hamiltonian 

if = ^(-<_iaX+i+Sa|), (5) 



where we assume periodic boundary conditions. The fact that 
crf_ia^af_^_i commute for different i and employing raising 

operator lJ, = cr^ — io'k-i'^k'^k+i allows to determine the en- 
tire spectrum of H for B = 0. The unique ground state of H 
for B = is the well-known cluster state [20,21], which has 
previously been studied as a resource in the context of quan- 
tum computation. It possesses a finite energy gap of AE = 2 
above its ground state [22]. For finite B the energy eigenval- 
ues of the system can still be found using the Jordan- Wigner 
transformation and a lengthy but straightforward calculation 
shows that the energy gap persists for ^ 1. The ex- 
act solution also shows that the system has critical points for 
\B\ = 1 at which the two-point correlation length and the en- 
tanglement length diverges. For any other value of B and in 
particular for B = the system does not exhibit a diverging 
two-point correlation length as is expected from the finite en- 
ergy gap above the ground state. Indeed, correlation functions 
such as 

^zz / 1 f^"" sinr . (L-l)r,s2 

Cff = (— / sm- — r-^dr) 

- , ^ + cos ^-^drr (6) 

can be computed and the corresponding correlation length 
can be explicitly determined analytically using standard tech- 
niques (see e.g. Fig. (2)) [23]. The two-point correlation 
functions such as Eq. (6) exhibit a power-law decay at the 
critical points = 1 while they decay exponentially for all 
other values of B in contrast to the anisotropic XF-model 
whose Cf£ correlation function tends to a finite constant in 
the limit of L — > oo for |S| < 1 [23]. This discrepancy is due 
to the finite energy gap the model in Eq. (5) exhibits above a 
non-degenerate ground state in the interval \B\ < 1. 

When we study three-spin interactions it is natural to con- 
sider the behaviour of higher-order correlations. For the 
ground state with magnetic field B = all three-point cor- 
relation except, obviously, {af_-^afaf_^-^) vanish. Indeed, if 
we consider n > 4 neighbouring sites and chose for each of 
these randomly one of the operators (Tx,(Jy, (Jz or 1 then the 
probability that the resulting correlation will be non- vanishing 
is given by p = 2~*^^+"^. For > however far more cor- 
relations are non- vanishing and the rate of non- vanishing cor- 
relations scales approximately as 0.858". This marked differ- 
ence which distinguishes B = is due to the higher symmetry 
that the Hamiltonian exhibits at that point. 

In the following we shall consider the localizable entangle- 
ment and the corresponding length as described in (ii). Com- 
pared to the two-point correlations, the computation of the lo- 
caUzable entanglement is considerably more involved due to 
the optimization process. Nevertheless, it is easy to show that 
the entanglement length diverges for B = 0. In that case 
the ground state of the Hamiltonian (5) is a cluster state with 
the property that any two spins can be made deterministically 
maximally entangled by measuring the cr^ operator on each 



3 



spin in between the target spins, while measuring the oper- 
ator on the remaining spins. Indeed, this property underUes its 
importance for quantum computation as it allows to propagate 
a quantum computation through the lattice via local measure- 
ments [20]. 

For finite values of B it is difficult to obtain the exact value 
of the localizable entanglement. Nevertheless, to establish a 
diverging entanglement length it is sufficient to provide lower 
bounds that can be obtained by prescribing specific measure- 
ment schemes. Indeed, for the ground state of (5) in the in- 
terval < 1 consider two spins 1 and i = 2fc + 1 where 
A: e N. Measure the operator on spin 2 and on all remain- 
ing spins, other than 1 and L, the operator. By knowing 
the analytic form of the ground state one can obtain the aver- 
age entanglement over all possible measurement outcomes in 
terms of the concurrence, that tends to Erx, = (l — |£?p)^''^ 
for k CO. This demonstrates that the localizable entan- 
glement length is infinite in the full interval \B\ < 1. This 
surprising critical behaviour for the whole interval \B\ < 1 is 
not revealed by the two-point or n-point correlation function 
which exhibit finite correlation lengths. For |_B| > 1 however, 
numerical results, employing a simulated annealing technique 
to find the optimal measurement for a chain of 16 spins, show 
that the localizable entanglement exhibits a finite length scale. 




Magnetic Field B 

FIG. 2. Both, the two-point correlation length for Ci£ (solid line) 
and the localizable entanglement length (dashed line) are shown for 
various magnetic field for chain of length 16. Note that the local- 
izable entanglement length diverges in the whole interval |_B| < 1 
while the two-point correlation length is finite in this interval. 

In Fig. (2) both the two-point correlation length and lo- 
calizable entanglement length are drawn versus the magnetic 
field. In the interval |_B| < 1 the entanglement length diverges 
while the correlation length remains finite. For finite tempera- 
tures the localizable entanglement becomes finite everywhere 
but, for temperatures that are much smaller than the gap above 
the ground state, it remains considerably larger than the clas- 
sical correlation length. This demonstrates the resilience of 
this phenomenon against thermal perturbations. 



To summarize, we have demonstrated that various Hamilto- 
nians describing three-spin interactions can be created in tri- 
angular optical lattices in a two-species Bose-Hubbard model. 
They can be realized in the laboratory with the near future 
cold atom technology. In fact, a study of the required experi- 
mental values reveals that with a tunneling coupling J/h ^10 
kHz [2] an experimentally achievable collisional coupling of 
U/h ~100 kHz is required. With these values a full nu- 
merical study demonstrates that the perturbative truncation is 
valid within a 4% error and a significant effect of the three 
spin interactions is obtained within the decoherence time of 
the system taken here to be 10ms. Previously, the systematic 
experimental creation of three-spin interaction Hamiltonians 
has been extremely difficult. The new capability for the sys- 
tematic creation of three-spin Hamiltonians and their possi- 
ble isolation from other interactions motivates the study of the 
properties of their ground states and here in particular of their 
phase transitions. Motivated by this we presented a particu- 
lar three-spin cluster Hamiltonian that exhibits a novel kind of 
critical behavior that is not revealed by two-point correlation 
functions. In addition, interactions such as erf (t|<7| presented 
here have proved to be of interest for quantum computation. 
They can implement multi-qubit gates, like the Toffoli gate, 
in essentially one step [24] reducing dramatically the experi- 
mental resources. 
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